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ABSTRACT 

We present a new method for fitting peculiar velocity models to complete flux 
limited magnitude-redshifts catalogues, using the luminosity function of the sources 
as a distance indicator. The method is characterised by its robustness. In particular, 
no assumptions are made concerning the spatial distribution of sources and their 
luminosity function. Moreover, selection effects in redshift are allowed. Furthermore 
the inclusion of additional observables correlated with the absolute magnitude - such 
as for example rotation velocity information as described by the TuUy-Fisher relation 
- is straightforward. 

As an illustration of the method, the predicted IRAS peculiar velocity model char- 
acterised by the density parameter [3 is tested on two samples. The application of our 
method to the TuUy-Fisher Marklll MAT sample leads to a value of /3 = 0.6 ± 0.125, 
fully consistent with the results obtained previously by the VELMOD and ITF meth- 
ods on similar datasets. Unlike these methods however, we make a very conservative 
use of the TuUy-Fisher information. Specifically, we require to assume neither the lin- 
earity of the TuUy-Fisher relation nor a gaussian distribution of its residuals. Moreover, 
the robustness of the method implies that no Malmquist corrections are required. 

A second application is carried out, using the fluxes of the IRAS 1.2 Jy sample 
as the distance indicator. In this case the effective depth of the volume in which the 
velocity model is compared to the data is almost twice the effective depth of the 
Marklll MAT sample. The results suggest that the predicted IRAS velocity model, 
while successful in reproducing locally the cosmic flow, fails to describe the kinematics 
on larger scales. 
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1 INTRODUCTION 

The study of the large-scale motions of galaxies in the Uni- 
verse may provide valuable information concerning the dy- 
namics of large-scale structures and the nature of the under- 
lying dark matter. According to the gravitational instability 
scenario, the peculiar velocity field (i.e. the deviation from 
the smooth Hubble flow) may be used to infer the power 
spectrum of the mass fluctuations on intermediate scales 
and to constrain the cosmological density parameter n (see 
for example Dekel 1994). 

Since the discovery of the Great Attractor (Lynden-Bell 
et al. 1988), the fleld has proven to be particularly active. 

Various observationnal programs have been completed, 
providing large and accurate datasets: e.g. the W91CL and 
W91PP samples (Willick 1990); MAT sample (Mathewson 
et al. 1992); HM sample (Han & Mould 1992); CF sam- 
ple (Courteau et al. 1993); Abell BOG sample (Lauer & 
Postman 1994); SCI sample (Giovanelh et al. 1997a); KLUN 



sample (Theureau et al. 1997); nearby SNIa sample (Riess 
et al. 1997); Marklll dataset (Willick et al. 1997b); SBF sur- 
vey (Tonry et al. 1997); SFI sample (Giovanelh et al. 1998); 
SMAC sample (Hudson et al. 1999); EFAR and ENEAR 
samples (CoUess et al. 1999 and Wegner et al. 1999); SCII 
sample (Dale et al. 1999); Shellflow survey (Courteau et al. 
1999); LPlOk survey (Willick 1999a). 

Another significant advance during the past decade has 
been the improved understanding of the statistical formal- 
ism underlying the use of galaxy distance indicators - and in 
particular the principles and practical methods of correct- 
ing for Malmquist bias. (See for example Hendry & Simmons 
1990, Teerikorpi 1990, BickneU 1992, Landy & Szalay 1992, 
Triay et al. 1994, WiUick 1994, Hendry & Simmons 1994, 
Sandage 1994, Willick et al. 1995, Freudling et al. 1995, 
WiUick et al. 1996, Rauzy & Triay 1996, Ekholm 1996, Triay 
et al. 1996, Rauzy 1997, WiUick et al. 1997b, GiovaneUi et 
al. 1997b, Theureau et al. 1998, Teerikorpi et al. 1999). 
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Several methods for extracting dynamical and kinemat- 
ical information from distance indicator datasets have been 
proposed, i.e. the POTENT method (Bertschinger & Dekel 
1989, Dekel et al. 1990, Bertschinger et al. 1990, Dekel et 
al. 1999 and references therein) and its variants (Rauzy et 
al. 1993 and 1995, Newsam et al. 1995); the ITF method 
(Nusser & Davis 1995, Davis et al. 1996, Da Costa et al. 
1998); the VELMOD method (Willick et al. 1997a, Willick 
& Strauss 1998). 

The comparison between the peculiar velocity or density 
fields inferred from distance indicator data with their corre- 
sponding fields derived from whole-sky redshift surveys has 
been one of the major issues addressed throughout the last 
decade. The question here is whether the spatial distribution 
of luminous matter e.g. the galaxies, traces the underlying 
mass fluctuations and if not, what are the properties of the 
"biasing" between the two fields? 

Up to now, the point has not received any consensual 
answer. Indeed, the application of POTENT to various dis- 
tance indicator datasets (Sigad et al. 1998 and references 
therein) favours a value of (3i — /bi ~ 1 for the linear 
"biasing" density parameter, while the VELMOD and ITF 
fitting methods lead to a value of /?/ ~ 0.5 (Davis et al. 
1996, Willick et al. 1997a, Riess et al. 1997, Da Costa et al. 
1998, Willick & Strauss 1998). The origins of this significant 
discrepancy have not yet been elucidated (see for example 
Strauss 1999, Willick 1999b). At least one of these methods 
is suffering from some systematic effects, i.e. some statisti- 
cal biases plaguing the estimate of the parameter /3 and not 
accounted for in the error analysis. 

This remark leads us to the object of the present paper. 
The POTENT, VELMOD and ITF methods all require, at 
some stage of the analysis, to assume some a priori working 
hypotheses concerning the characteristics of the distance in- 
dicator dataset. For TuUy-Fisher data for example, it will 
be assumed that the TuUy-Fisher law is well described by a 
linear relation. These methods apply moreover under the hy- 
pothesis that the observational selection effects obey partic- 
ular conditions. How are the results affected if one or more of 
the working assumptions fails to be satisfied by the dataset 
is generally a question not addressed in the error analysis. 

The philosophy of the method we present herein is to 
reduce as far as possible the number of a priori hypotheses 
concerning the distance indicator sample. A direct conse- 
quence is that the range of application of the method will 
be considerably broadened. 

The statistical background of the method is presented 
section 2. Its potential is illustrated by testing the predicted 
IRAS peculiar velocity model on two samples. In section 3, 
we perform the analysis using the fluxes of the IRAS 1.2 
Jy survey as the distance indicator. We have deliberately 
chosen this sample looking to demonstrate the wide range of 
application of the method. Where the POTENT, VELMOD 
or ITF methods would not have been successful in extracting 
kinematical information from this dataset, our method does. 
We also treat a more classical case, the TuUy-Fisher Marklll 
MAT sample, in section 4. Finally, in section 5 we summarise 
our conclusions. 



2 THE METHOD 



2.1 The statistical model 

It is assumed herein that the distribution function of the 
absolute magnitudes M of the population, i.e. the luminosity 
function f{M), does not depend on the spatial position r = 
{r,l,b) of the galaxies. The probability density describing 
the sample splits in this case as 

dP oc dPr X dPM = p{r, I, b) cosbdldbdr x f{M)dM (1) 

where p{r,l,b) is the spatial distribution function of the 
sources. 

The application of the method will be restricted to 
samples strictly complete up to a given magnitude limit 
miini, or in other words where the selection function in ap- 
parent magnitude is well described by a sharp cut-off, i.e. 
tp{m) = 0(miini — m) with 6{x) the Heaviside function. Ac- 
counting for selection effects, the probability density of the 
sample may be rewritten as 

dP= ^h{fj.,l,b)cosbdldbdfif{A4)dMe{mu,-n-m) (2) 

where /i = m — M = 5 log^ r -f- 25 is the distance modulus, 
h{fi, I, b) the line-of-sight distribution of ^ and A is the nor- 
malisation factor warranting J dP = 1 . For convenience in 
notation the angular dependence in I and b will be hereafter 
implicit. Observational selection effects in apparent magni- 
tude then introduce a correlation between M and fi. 

The milestone of the method consists in defining the 
random variable ( as follows 



F{M) 



(3) 



where F{M) = j^' f{x)dx stands for the cumulative dis- 
tribution function in M and = Mnj^ifJ') is the maxi- 
mum absolute magnitude for which a galaxy at distance ^ 
would be visible in the sample (e.g. Miiin(/x) = miim — fi if 
the k-correction is neglected). The volume element may be 
rewritten as 



(4) 



and by definition the random variable ^ for a sampled galaxy 
belongs to the interval [0, 1]. The probability density of Eq. 
(^) reduces thus to 



dP = jhin) F(Mnm{fi)) d^i X 6»(C) 6»(1 - C) 



(5) 



with A — J h{fi) F{Mii^{fi)) dfi. Note that the probability 
density dP^^ — ^h{fj,) _F(Miim(A')) dfj, describes the observed 
spatial distribution function of the sources. It follows from 
Eq. (||) that: 

• PI: is uniformly distributed between and 1. 

• P2: ( and fi are statistically independent, i.e. the dis- 
tribution of does not depend on the spatial position of the 
galaxies. 

Property PI may be used to construct a test for assessing 
the completeness of the sample in apparent magnitude. The 
details of this test are presented in a separate paper (Rauzy, 
in preparation), although we apply it to the Mark III MAT 
sample later in this paper. The new method we propose 
hereafter for fitting peculiar velocity field models is based 
on property P2. 
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2.2 Estimate of the random variable 

The random variable ^ can be estimated without any prior 
knowledge of the cumulative luminosity function F{M). To 
each data point with coordinates [Mi,^i) is associated the 
region Si = SiVJ S2 defined as 

• Si = { (M, ^l) such that M < Mi and /J < /li } 

• 5*2 = { (M, ^l) such that Mi < M < Ml^^-^-^ and a* < M» } 

The random variables M and p are independent in each sub- 
sample Si since by construction the cut-off in apparent mag- 
nitude is superseded by the constraints M < Ml^^ and /i < 
/ii (see figure 1). This implies that the number of points 
belonging to Si is proportional to fi^) dM — F{Mi), 
the number of points Ui in Si = Si U S2 is proportional to 
F{Mii^) and that the quantity 



= 



rii + l 



(6) 



is an unbiased estimate of the random variable Equiv- 
alently the estimator d may be defined as the normalised 
rank of the point Ali when the M's are sorted by increasing 
order within the subsample Si (see Efron & Petrosian 1992) . 



2.3 Radial peculiar velocity field models 

Let us first assume that the true radial peculiar velocity field 
v{r) can be described by a one-parameter velocity model 
vp{r), i.e. there exists a solution (3* satisfying v/s* (r) = v{r). 

For a given value of the parameter /?, the model de- 
pendent variables fijs and M/s can be computed (modulo the 
value of the Hubble constant Hq) from the observed redshift 
z and apparent magnitude m following 
cz 

^1/3= 5 logio — + 25 -Up ■ Mp = m- fif} (7) 

-no 



where the quantity up is defined as 



Uf3 



-5 log 



(8) 



The quantities fip and Mp are related to the true absolute 
magnitude M and distance modulus fi via 



^1/3=^1 + Up* 



Up 



Mp ^ M - up* +up 



(9) 



Computing (p from fip and Mp as 
for the probability density of Eq. (H), 

dP = jh{fi)FiMii^{fip))dfi Cp e{Cp)eii 



Droposed in Eq. (H) gives 



C/3)rfC/3 (10) 



where Cp takes the following form when {up* — up) <C 1 (or 
equivalently {vp* — vp) <^ cz), 



f{M) 
fiMp) 



l + {up-up*) (In f)'{Mp) 



(11) 



Because the absolute magnitude Mp, and hence the quantity 
(In/) (Mp), is correlated with the random variable (^p, Cp 
acts as a correlation coefficient between (^p and the proposed 
velocity field model up when j3 ^ l3* . On the other hand if 
j3 — P* these quantities are statistically independent, since 
in this case, according to property P2, (p = does not 
depend on the spatial position of galaxies and therefore on 
any function up{r). 



It thus turns out that any statistical test of indepen- 
dence between (^p and up provides us with an unbiased es- 
timate of the value of /3*. In particular the coefficient of 
correlation p{C,p,up) has to vanish when (3 = i.e. 



/3 = /3* 



p{C,p,up) = 



(12) 



As revealed by Eq. (|ll]), the accuracy of this estima- 
tor is related to the amplitude of the correlation between 
{\a f)' {Mp) and ^^3. The steeper the function (In/)', or in 
other words the smaller the dispersion of the luminosity 
function f{M), the more accurate is the estimate of the ve- 
locity parameter l3, as expected. In practice, this accuracy 
can be obtained through numerical simulations by analysing 
the influence of sampling fluctuations on the coefficient of 
correlation p{C,p,up). 

The presence of a small-scale velocity dispersion (say 
of amplitude ct„), not described by the velocity model vp, 
introduces according to Eq. (^) a correlation between the 
derived quantities fj,p and Mp, and consequently between 
the variables pp and C,p. Anyway since it is the correlation 
between the velocity model up and C,p which is considered 
herein, and because the random velocity noise is not sup- 
posed to be correlated with up, the presence of a small-scale 
velocity dispersion is not expected to drastically bias the es- 
timator proposed Eq. (|l^), at least as long as the variations 
of the quantity up{r) axe smooth at the scale cr„. 

Thanks to the introduction of the random variable C,, 
an unbiased estimate of the parameter (5 has indeed been 
obtained using a null-correlation technique. Null-correlation 
approaches are characterised, in general, by their robust- 
ness - i.e. some of the functions entering the statistical 
model are not required to be fully specifled (see for exam- 
ple Fliche & Souriau 1979, Bigot et al. 1991, Triay et al. 
1994, Rauzy 1997). Unlike the maximum likelihood meth- 
ods, e.g. the method proposed by Choloniewski (1995) and 
the VELMOD method of Willick et al. (1997a), no a priori 
assumptions have been made here concerning the speciflc 
shape of the luminosity function and the spatial distribution 
of the sources. In particular homogeneous as well as inhomo- 
geneous Malmquist biases are automatically accounted for. 
Note also that selection effects in distance or redshift are 
allowed since Eq. (^^ accepts any extra terms of the form 

Tp{il,,Up). 



2. 3. 1 Orthonormal decomposition of the velocity field 

It is worthwhile to mention that, with reg ard to Eq. (0), the 
orthonormal decomposition of the velocity fleld, proposed in 
the ITF method (Nusser & Davis 1995) , may be also applied 
herein. To see this we proceed as follows. 

It is assumed hereafter that there exists a N-dimensional 
vector /3* = {(3i, /^l > ••■> Pn) such that the quantity u(r) can 
be decomposed as 



u{r) = up*{r) 



N 



(13) 



where ui{r), U2{r), itjv(r) is a set of functions verifying 
the following orthonormality condition. 



Cov(ui, Uj) 



(14) 
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Figure 1. Distance modulus versus absolute magnitude for the 
60 /^m IRAS 1.2 Jy sample (5321 galaxies). The procedure for 
estimating the random variable C, from Eq. (6) is illustrated. 



with 5ij the Kronecker symbol and Coy {ui,Uj) the covari- 
ance of Ui and Uj on the sample. Note that within the ap- 
proximation P*Uj{r) <^ 1 (or equivalently PjVj{r) <^ cz), 
Eq. (^) implies that 



.(r)~^/3>,(r) 



(15) 



The coefficient C^g introduced in Eq. ( [11^ rewrites, as long 
as [Pj — P*)uj{r) -C 1 is satisfied, as 



~ 1 + (In/)' {Mp) X ^(/3, -/3*)^,(r) 



(16) 



o Y 

C 
o 

X2 O 
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which implies that for each function Ui{r) 



-24 -22 -20 -18 -16 

Absolute magnitude M 

Figure 2. Reconstruction of the Cumulative Luminosity Func- 
tion of the IRAS 1.2 Jy sources using the C~ method of Lynden- 
Bell (1971). The reconstructed CLE depends on the assumed ve- 
locity model. 



p{u., Cp) cc Y^{f3j - /?;) Cov(m„ Uj) (17) 

and thus, because of the orthonormality condition of Eq. 
(0), that 



A = A* 



p{Cl3,Ui) = 



(18) 



which provides us with the statistically independent esti- 
mates of the A'^ parameters PI, (52, Pn ■ The procedure 
for constructing the orthonormal family iii(r), U2{v), 
Mjv(r) from an arbitrary set of A'^ independent functions is 
described in Nusser & Davis (1995). 



3 APPLICATION TO THE IRAS SAMPLE 

The method described above is herein applied to the 60 
Hm IRAS 1.2 Jy sample (Fisher et al. 1995), a magnitude- 
redshift catalogue containing 5321 galaxies and complete up 
to a flux Feo = 1.2 Jy. Distance modulus and absolute mag- 
nitude are computed using Ho = 100 km s"^ Mpc"^ for the 



value of the Hubble constant (the cut-off in apparent magni- 
tude is expressed as mum = 14.3187 with this notation). The 
magnitudes have been k-corrected assuming a spectral slope 
a — —2, implyingthat the maximum absolute magnitude in- 
troduced in Eq. ij^ reads as MnmC/i) = "^lim — kcor(/^) — M. 
The distribution of the sources in the M-fj, plane is shown 
figure 1. 

The peculiar velocity field model tested is the one- 
parameter predicted IRAS velocity field characterised by the 
parameter P = Q,o^/bi (Strauss et al. 1992). 



3.1 Using apparent magnitude as the distance 
indicator 

We have first applied the method using the apparent magni- 
tude of IRAS galaxies as the distance indicator (restricting 
the analysis to the 4115 galaxies within the redshift range 
1000-12000 km s~^). The Cumulative Luminosity Function 
(CLF) F{M) is presented in figure 2. Here the CLE has been 
reconstructed using the C~ method of Lynden-Bell (1971) 
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(^=0.6) 

Figure 3. Scatterplot indicating the correlation between tlie ve- 
locity modulus Ufj and the random variable (^p for /3 = 0.6. A 
observed sample correlation coefficient of p(C/3i"/3) — 0.02 was 
found in this case. 



for two different velocity models, (3 — corresponding to 
the case where peculiar velocities are neglected and /3 — 1. 
Note that the presence of peculiar velocities does not affect 
drastically the CLF reconstruction. The luminosity function 
of the IRAS galaxies does not exhibit any turnover towards 
the faint-end tail, at least within the observed range of mag- 
nitudes. This means that, because of the large dispersion of 
such a distance indicator and even if the number of galaxies 
is large, one cannot expect very strong constraints on the 
velocity model tested. 

The correlation between the random variable C^p and 
the velocity modulus up for (3 = 0.6 is shown in figure 3. 
Variations of the coefficient of correlation as a function of the 
parameter /3 are given in figure 4. This curve is a monotonic 
function, as expected. The preferred value of /3 is the one 
corresponding to p{(,i3,up) = (here /3 = 0.1 — 0.15). 

Monte Carlo simulations have been used to calculate 
the discrepancy between p{(i3,up) and zero due to sam- 
pling fluctuations. Each simulation is a sample containing 
-^gai = 4115 galaxies with the same fip and up as observed 
and for which the random variable (^p is computed following 
Eq. (^) where the rank is randomly generated according 
to a discrete uniform distribution between 1 and rn. The cu- 
mulative distribution function of p{Cf3,up) obtained from a 
large number of simulations, under the null hypothesis that 
the true value of p{(p,up) — 0, is shown in figure 5. Note 
that in practice this distribution does not depend on the 
amplitude of the quantity up since the coefficient of correla- 
tion is a scale-free estimator of the correlation between two 
random variables. 

The cumulative distribution function allows us to eval- 
uate the probability that the observed p is by chance greater 
than a given value, due to sampling fiuctuations. A one-sided 
rejection test for the f3 parameter can thus be constructed. 
For example, the models with P > 0.7 can be rejected with 
a confidence level of 95% and /3 > 1.1 with a confidence 
level of 99%. So our method applied to the IRAS sample, 
using the apparent magnitude of the galaxies as a distance 
indicator, permits already at this stage to reject high values 
for the parameter /3. 
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Figure 4. The observed coefficient of correlation p(C/3,«/3) as a 
function of the parameter Confidence levels of rejection for the 
parameter /3 are calculated from figure 5. 
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Figure 5. Cumulative distribution function of p due to sampling 
fiuctuations. This curve has been obtained by Monte Carlo sim- 
ulations (see text). 



3.2 Introduction of a second parameter 

In a second step, the analysis is refined by taking into 
account the observed correlation between the absolute 
magnitude M and some "colour index" defined as p = 
2.51ogj^Q(Fioo/f6o) (with Fioo the fiux at 100 /im). The data 
have been grouped in 8 classes by interval of p (see figure 
6). Because of the (weak) correlation between p and M, the 
spread of the luminosity function for each of these classes 
taken individually is expected to be slightly smaller than 
the spread of the global luminosity function, and thus the 
accuracy of the distance indicator somewhat improved. Fig- 
ure 7 illustrates such a trend. 

For each individual class, the random variable C,p is 
computed according to Eq. (^). The correlation between C,p 
and the velocity modulus up is after that evaluated for the 
whole sample. The influence of sampling fluctuations is es- 
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Figure 6. Observed correlation between tiie "colour index" p = 
2.51ogio(-Pioo/.F6o) and the absolute magnitude M. The sample 
has been subdivided into 8 classes, C, according to the value of 
P- 
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Figure 7. reconstruction of the Cumulative Luminosity 

Function for the 8 individual d classes in p (dark curves). The 
global CLF (i.e. no class) of figure 2 is shown for coriiijarison (grey 
curve). The CLF's have been arbitrarily normalised to 10~* at 
M = -20. 



timated from Monte Caxlo simulations, as described in the 
previous section. The results are presented in figure 8 in 

terms of the confidence level of rejection for the parame- 
ter 13. More explicitly, the quantity plotted in ordinate is 
1 — 2Prob(p < — |/5obs(/3)|) whore the probability that the 
coefficient of correlation p is less than — |pobB(/3)| due to sam- 
pling fiuctuations is given by the cumulative distribution 
function of p. 

The method was first applied to the galaxies in the 
redshift range 1000-12000 km s"^ It is found that /? € 
[—0.35, 0.25] at Icr, and that models with /3 > 0.5 can be 
rejected with a confidence level of 95%. This result is in 
disagreement with most of the analyses based on TuUy- 
Fisher data e.g. VELMOD on Marklll (Wilhck & Strauss 
1998), ITF method on SFI (Da Costa et al. 1998), ROBUST 
method on Marklll MAT sample (see next section), favour- 
ing a value of /3 ~ 0.5. We interpret this discrepancy as 
follows. 

When fitting a velocity model to data, the natural 
weight assigned by the fitting procedure to each galaxy is 
roughly proportional to the inverse of its redshift, because 
the accuracy of the distance indicator decreases as 1/z. The 
mean effective depth of the volume where the velocity model 
is compared to data has to be estimated using these weights. 
For our first sample with z € [1000, 12000] km s~^, we find 
a mean effective depth of 3800 km s~^ (see figure 9). 

In order to mimic the effective volume sampled by 
TuUy-Fisher data wo have applied the method to a trun- 
cated version of the IRAS sample containing 1621 galaxies 
with z € [500, 5000] and galactic latitude \b\ > 20 (the mean 
effective depth of this sub-sample is now 2200 km s~^, sec 
figure 10). Figure 8 shows that the value of (3 estimated 
from this truncated sample is fully consistent with the val- 
ues obtained using TuUy-Fisher data. An interpretation of 
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Figure 8. Confidence level of rejection for the parameter (3 (see 
text). 



these results could bo that the predicted IRAS velocity field 
model, while successful in reproducing locally the cosmic 
fiow, fails to describe the kinematics on larger scales. 

However, as pointed out by our anonymous referee, the 
results derived above apply only to the extent that the pho- 
tometry of the IRAS sample does not suffer from systematic 
errors. It is worthwhile to stress again that the philosophy of 
the ROBUST method is to impose very few assumptions - 
only that the luminosity function is independent of position, 
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Weighted: <z> = 3800 km/s 
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Figure 9. Redshift distribution of the IRAS galaxies with z G 
[1000, 12000] km s"^ (4115 galaxies). The dark histo gram rep- 
resents this distribution when accounting for the natural weight 
oc 1/z. 



Weighted: <z> = 2200 km/s 
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Figure 10. Redshift distribution of the IRAS galaxies with z G 
[500, 5000] km s"! and galactic latitude |6| > 20 (1621 galaxies). 
The grey histogram accounts for the natural weight oc l/z. 
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Figure 11. The TuUy-Fisher relation for the Marklll MAT sam- 
ple (Willick et al. 1997b). A subsample of 318 galaxies (dark dots) 
complete in apparent magnitude up to muni = 11.25 mag. has 
been extracted. This subsample is divided in 6 classes Cj con- 
taining approximately 50 galaxies each. 



goal here is to present a now method, we feel that such an er- 
ror analysis is unwarranted. It is also worth noting that sys- 
tematic errors in the IRAS photometry - if indeed present 
- would also affect methods for reconstructing the IRAS 
predicted peculiar velocity model. In particular, they would 
generate some systematic variations in the spatial selection 
function entering the weighting scheme used in density and 
velocity reconstruction (Strauss et al. 1992, Branchini et al. 
1999), leading to systematic discrepancies between the IRAS 
predicted velocity field and the true cosmic flow. Moreover, 
it is interesting to note that the ROBUST method could 
easily be extended to provide a non-parametric test of the 
uniformity of photometry in redshift surveys; we will con- 
sider such an application in future work. 



that the sample is strictly complete in apparent magnitude 
and that redshift and apparent magnitude measurements 
are not affected by systematic biases. Thus, our analysis 
could be affected if the IRAS photometry were (mildly) non- 
uniform. To evaluate the amplitude of such effects on our re- 
sults would essentially require to adopt a realistic model for 
these systematic photometry variations. Since our primary 



4 APPLICATION TO THE MARKIII MAT 
SAMPLE 

We have shown in the previous section that our method per- 
mits to extract valuable information on the cosmic velocity 
field by using the fluxes of the IRAS galaxies as a distance 
indicator. The potential of our method is illustrated in this 
section by treating a more classical case: the TuUy-Fisher 
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Figure 12. For each of the 6 classes d, the reconstructed 
Cumulative Luiiiiiiosity Function (dark dots). The solid curves 
are the expected cumulative distribution function of M assuming 
a linear Direct (i.e. Forward) TuUy-Fisher relation with Gaus- 
sian residuals and calibrated using the values proposed in Willick 
et al. (1998) (see text). For convenience in comparison, the C~ 
reconstructed CLF's Frec(JW") have been normalised such that 
Frec(A/i) = 0.5 at Mi the mean absolute magnitude for each 
classes Cj (dotted lines). 
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Figure 13. Confidence level of rejection for the parameter /? 
(/3 = 0.6 ±0.125). 
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Figure 14. Redshift distribution of the 318 Marklll MAT galax- 
ies with z > 500 km s~^ and apparent Tband magnitude m < 
11.25 mag. The dark histogram takes into account the natural 
weight oc 1/2. 



Marklll MAT sample (Willick et al. 1997b and references 
therein) . 

In a first stop wo have scloctod from the 1355 galaxies 
of the Marklll MAT catalog a subsample for which the se- 
lection effects in apparent magnitude are well described by 
a Heavisido cut-off, 9(m — mum). Assuming that the galaxies 
are homogeneously distributed in space, the value of mum 
can be found by analysing the variations in the logarithm of 
the cumulative count as a function of the limit in apparent 
magnitude - see for example figure 4 in Rauzy (1997). In 
a separate paper (Rauzy, in preparation) we propose any- 
way a simple tool for assessing the completeness in apparent 
magnitude of redshift-magnitude catalogues which docs not 
require any cissumption concerning the spatial distribution 
of the sources. This method is closely based on the statis- 
tical test presented in Efron & Petrosian (1992). We have 
apphed this test to the Marklll MAT sample, finding that 
the completeness in I-band corrected apparent magnitude is 
satisfied up to mnm = 11.25 magnitudes. 

Discarding, in addition, the galaxies with z < 500 km 
s~^, we are left with a subsample containing 318 galaxies. 
Note that the selection effects in redshift which affect the 
Marklll MAT sample (Willick et al. 1996) are not a source 
of problems for our method. The Tully-Fisher relation for 
the 318 extracted galaxies is shown in figure 11. 

In a second step, we divided the subsample in 6 classes 
C'i according to the value of r/, from the slower rotators in Ci 
to the faster rotators in Ce (see figure 11). The C~ recon- 
structed cumulative luminosity function is shown for each 
class in figure 12. If the binning in the log line-width pa- 
rameter ri were narrow enough, one would expect the recon- 
structed luminosity function for the class Ci to be centered 
on Mi = aoTFTji + 6dtf and of dispersion ctdtf (»?<), where 
rji =< V >Ci denotes the mean of the rj's in d and odtf, 
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feoTF and CTDTF('7i) are respectively the slope, zero-point and 
dispersion of the Direct (i.e. Forward) TuUy-Fisher relation. 
We plot in figure 12 the expected culmulative luminosity 
function when a linear DTF relation with gaussian residu- 
als is assumed and using the values proposed in Willick et 
al. (1998) for the calibration parameters. The width of the 
bins has been accounted for by adding in quadrature to the 
dispersion aoTFiVi) the product of IodtfI and the standard 
deviation of the rj's in each class d . For convenience in 
comparison, Each reconstructed CLF Frec{M) has been 
normalised such that Frec(Afi) = 0.5. 

It is not clear from figure 12 whether a linear DTF rela- 
tion with gaussian residuals is plainly successful in reproduc- 
ing the data. Note in particular that the reconstructed 
CLF's for the slow rotators do not exhibit a turnover to- 
wards the faint magnitudes. To question the validity of a 
linear TuUy-Fisher relation is anyway beyond the scope of 
the present paper. At this stage it is however worthwhile to 
mention that our method makes a very conservative use of 
TuUy-Fisher information. In fact, we require to assume nei- 
ther a linear TF relation nor a gaussian distribution for its 
residuals. 

The result of the ROBUST method applied to the 
Marklll MAT catalog is shown in figure 13. The analysis 
was performed as described in section 3.2. The mean ef- 
fective depth of the subsample is 2100 km s~^ (see figure 
14). We find a value of /3 = 0.6 ± 0.125, in complete agree- 
ment with the VELMOD and ITF method applied to similar 
Tully-Fisher data (WiUick et al. 1997a and 1998, Davis et 
al. 1996, Da Costa et al. 1998). 

We want however to emphasise the robustness of our 
approach compared to these two fitting methods. Firstly, no 
assumptions have been made herein concerning the linearity 
of the Tully-Fisher law, as is required by both the VELMOD 
and the ITF method. Secondly, we do not need the sample 
to be free of selection effects in the log line- width parameter 
rj as is the case for the ITF method. Thirdly, the spatial 
distribution of the sources, the selection effects in redshift 
and the shape of the distribution function of the TF resid- 
uals need not be specified, as is required by the maximum 
likeUhood VELMOD method. 



5 CONCLUSION 

We presented a method for fitting peculiar velocity mod- 
els to complete fiux limited magnitude-redshift catalogues, 
using the luminosity function of the sources as a distance 
indicator, i.e. assuming that the distribution function of the 
absolute magnitudes of the galaxies does not depend on the 
spatial position. 

Our method is based on a null-correlation approach. 
For a given peculiar velocity field model parametrised by a 
parameter /?, we defined a random variable (p, computable 
from the observed redshifts and apparent magnitudes of the 
sampled galaxies, which has the property of being statis- 
tically independent on the position in space (and thus on 
the modelled radial peculiar velocities themselves) if and 
only if the parameter ,3 matches its true value /3*. Therefore 
any test of independence between the random variable C/3 
and the modelled velocities or similar quantities provides us 
with an unbiased estimate of the value of (3*. The method 



can be easily generalised to velocity models parametrised by 
an A''-dimensional vector /3 = {/3i,f32, ...,/3iv). 

The method is characterised by its robustness. No as- 
sumptions are made concerning the spatial distribution of 
sources and their luminosity function and selection effects in 
redshifts are also allowed. The required strict completeness 
in apparent magnitude can moreover be checked indepen- 
dently (Rauzy, in preparation). Furthermore the inclusion 
of additional observables correlated with the absolute mag- 
nitude is straightforward. 

The predicted IRAS peculiar velocity model charac- 
terised by the density parameter /3 has been tested on two 
samples, the TuUy-Fisher Marklll MAT sample and the 60 
/im IRAS 1.2 Jy sample using the fluxes as the distance 
indicator. 

The application of our method to the Marklll MAT 
sample gives a value of /3 = 0.6 ± 0.125, in excellent agree- 
ment with the results obtained previously by the VELMOD 
and ITF methods on similar datasets. Our method is how- 
ever more robust than these two fitting methods. In partic- 
ular, we make a very conservative use of the Tully-Fisher 
information. We do not require to assume the linearity of 
the Tully-Fisher relation nor a gaussian distribution of its 
residuals. 

We showed that our method allows to extract some 
valuable informations on the peculiar velocity field from the 
fiuxes of the IRAS 1.2 Jy sample. The poor accuracy of the 
distance indicator (due to the broad spread of the luminosity 
function) is balanced in this case thanks to the large num- 
ber of galaxies contained in the sample. The IRAS sample 
permits to probe the cosmic fiow at larger scales. Indeed, 
the mean effective depth of the volume in which the veloc- 
ity model is compared to the data is almost twice the mean 
effective depth of the Marklll MAT sample. 

The application of our method to an IRAS subsample 
truncated in distance, of an effective depth similar to the 
Marklll MAT sample, gives a value of /3 in accord with the 
values obtained using Tully-Fisher data. On the other hand 
when the application is performed on the whole sample, we 
found that the predicted IRAS velocity models with /3 > 0.5 
can be rejected with a confidence level of 95%. These re- 
sults suggest that the predicted IRAS velocity model, while 
successful in reproducing locally the cosmic flow, fails to de- 
scribe the kinematics on larger scales. 

Note that these results do not lead to dismiss the linear 
■'biasing" paradigm. As the errors on the predicted IRAS 
velocity fleld increase with distances, it could be that the 
predictions at the scales considered herein, i.e. beyond 5000 
km s~^, drastically differ from the true cosmic flow (see for 
example Davis et al. 1995). 
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